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Abstract 

The QCD one-loop renormalization is restudied in a mass-dependent 
subtraction scheme in which the quark mass is not set to vanish and 
the renormalization point is chosen to be an arbitrary time-like mo- 
mentum. The correctness of the subtraction is ensured by the Ward 
identities which are respected in all the processes of subtraction. By 
considering the mass effect, the effective coupling constant and the 
effective quark mass are given in improved expressions which are dif- 
ferent from the previous results. 
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1. Introduction 

The asymptotically free property of Quantum Chromodynamics (QCD) has 
been widely applied to analyze the scaling behaviors of high energy processes 
and attracted much attention to investigate the QCD renormalization^"^^. 
This property was at first discovered in Refs.(l) and (2) from studying the 
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one-loop renormalization and the solution of renormalization group equa- 
tion (RGE)^''~^^satisfied by the coupling constant. The anomalous dimen- 
sion ( or the /3-function) appearing in the equation was computed from the 
renormalization constants which are ordinarily given by mass-independent 
subtractions such as the minimal subtraction (MS)®, the modified minimal 
subtraction''and the momentum space subtraction(MOM)^ . However, there 
appears a serious ambiguity problem^^ that different subtraction schemes 
give different results for a renormalized quantity in a finite order perturbative 
calculation. In the MS schemes which is carried out in the dimensional reg- 
ularization procedure, only the divergent terms proportional to the factor - 
(where e = 2 — ^) are subtracted. The MS scheme was demonstrated to give 
the perturbative expansion in the coupling constant for a physical quantity 
which is of worse convergence. To improve the convergence, the MS scheme 
was proposed to subtract the divergent part of a Feynman diagram which 
contains the factor k = ^ + In 47r — 7 where the unphysical terms In 47r — 7 
arise from a special way of analytical continuation of the dimension from n 
to 4^. The MS and MS schemes are simple and respect the Ward identity, 
but considered to be nonphysical^. In the MOM scheme, the divergence in 
a Feynman integral is subtracted at an Euclidean momentum, pf = —/i^. 
This subtraction was viewed as physical and shown to give the perturbative 
series of better convergence. Nevertheless, the renormalization constants ob- 
tained in this scheme do not automatically satisfy the Ward identity. To solve 
the ambiguity problem, several prescriptions were proposed in the past^°~^^. 
In our preceding paper on the QED renormafization^^, we treated the am- 
biguity problem from a different angle. It was argued that a subtraction 
scheme must respect necessary physical and mathematical principles shcu as 
the gauge-symmetry (the Ward identity), the Lorentz-invariance( the energy- 
momentum conservation) and the mathematical convergence. In addition, it 
was proved that the solution of the renormalization group equation (RGE) 
for a renormalized quantity (for example, a wave function, a propagator or a 
vertex) can be uniquely determined by applying the renormalization boudary 
condition without any ambiguity. As a result, an exact S-matrix element can 
be expressed in the form as given in the tree-diagram approximation ex- 
cept that the coupling constant, the fermion mass and the gauge parameter 
become effective (running )ones. Therefore the tesk of renormalization is 
reduced to find these effective quantities. 

In the mass-independent subtractions mentioned before, the fermion mass 
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is set to be zero, or say, only the massless fermion propagator ^ is used in 
the calculation of divergent Feynman diagrams. The rationality of the mass- 
independent subtraction was argued as foUows^'^^. In the conventional per- 
turbation series which is expanded in coupling constant, the massive fermion 
propagator can be expanded as such a series 



Substitution of this series into the coventional perturbation series leads to a 
new perturbation series in which the fermion propagator becomes massless 
and the fermion mass, like the coupling constant, can be treated as another 
expansion parameter. Working with this kind of perturbation series, one may 
perform the mass-independent subtraction. It is obvious that in order to ob- 
tain a quantity which is of a certain order of coupling constant in the sense of 
the conventional perturbation theory, according to Eq.(l.l), one has to com- 
pute an infinite set of perturbative terms in the new perturbation theory. In 
the large momentum limit, the massive propagator is reduced to the massless 
one. In this case, there will be no difference between the both perturbation 
theories mentioned above. However, for a process whose energy is not too 
high, particularly, in the energy region near the meson production threshold, 
the mass in the fermion propagator can not be negligible. In this case, to see 
the effect of fermion mass on the renormalization , it is convenient to carry 
out a mass-dependent renormalization in the framework of the conventional 
perturbation theory as was done in the previous literature^'^''. Obviously, 
the MOM scheme is suitable for this kind of renormalization. In Ref.(3), the 
subtraction was carried out at an arbitrary Euclidean point, p'f = —fi'^ and 
in the Landau gauge, while, in Ref.(lO), the renormalization was done at an 
off mass shell point, p} — m} = —jj?' and in the Feynman gauge. Both of the 
subtractions give the same one-loop anomalous dimension of the coupling 
constant, but different anomalous dimensions of the quark mass. 

In this paper, we wish to restudy the QCD mass-dependent renormaliza- 
tion in the MOM scheme along the line as described in our preceding paper. 
The new features of this study are: (1) The subtraction exactly respects 
necessary physical and mathematical principles. Owing to the restriction of 
these principles, the renormalization will be put on the faithful basis and 
have no ambiguity. For example, by the convergence principle, we arc not 
allowed to employ the divergent form of renormalization constants to do a 
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meaningful calculation because a divergent quantity is not well-defined math- 
ematically and therefore can not be unambiguously evaluated by any compu- 
tational rule. The correct procedure of computing the anomalous dimension 
in a RGE is usage of regularized froms of the renormalization constants. 
By this procedure, a renormalized quantity given by solving its RGE will 
be uniquely determined by its anomalous dimension and boundary condi- 
tion; (2) the renormalization point is taken to be a time-like (Minkowski) 
momentum. The subtraction with this renormalization point will be called 
generalized mass-shell scheme(GMS) because this scheme can naturally lead 
to the results given in the mass-shell subtraction scheme. It will be shown 
that the behavior of a renormalized quantity derived in the GMS scheme is 
different from the ones obtained in other subtraction schemes; (3) The QCD 
one-loop effective coupling constant and effective quark masses derived in the 
GMS scheme and in general gauges are given rigorous and explicit expres- 
sions. These expressions will go over to the results given in the MS scheme 
in the large momentum limit. 

The remainder of this paper is arranged as follows. In Sect. 2, we will 
describe the derivation of the renormalization constant of coupling constant. 
This renormalization constant is chosen to be determined by the renormal- 
ization constants of the gluon propagator, the ghost particle propagator and 
the ghost vertices. The renormalization constant of gluon propagator is given 
based on the Ward identity. It will be shown that the other two renormaliza- 
tion constants for the ghost particle propagator and the ghost vertex will give 
their anomalous dimensions which are renormalization-scheme-independent. 
In Sect. 3, we will derive an explicit expression of the effective coupling con- 
stant by solving its RGE and show its asymptotic behavior. In Sect. 4, we 
will give a derivation of the Ward identity satisfied by the quark-gluon vertex 
and prove that in the approximation of order g^, this identity is reduced to 
the form as found in QED and thus the subtraction version of the quark self- 
energy is similar to that for the electron self-energy. In Sect. 5, an explicit 
expression of effective quark masses will be given by solving its RGE. We 
will end this paper with some comments and discussions. 
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2. The Renomalization Constant of Coupling 
Constant 



The renormalization constant of coupling constant Zg is defined by 

9 = Zggn (2. 1) 

where g^i and g denote the rcnormahzcd and unrenormahzed couphng con- 
stants respectively. According to the Ward identity'^°, the Zg can be ex- 
pressed in different ways. In this section, we prefer to use the following 
expression 

= ^ (2. 2) 

where Z3, Z3 and Z\ are the renormalization constants for the gluon prop- 
agator, the ghost particle propagator and ghost vertex respectively. In the 
following, we will describe how they are determined by the subtractions of 
the gluon self-energy, the ghost particle self-energy and the ghost vertex cor- 
rection in the GMS scheme. The one-loop diagrams for the gluon self-energy, 
the ghost particle self-energy and the ghost vertex correction have already 
been calculated in the literature^'^^'^^ by the dimensional regularization . So, 
we may directly quote the results and put some emphases on essential points 
of the subtraction procedure in the GMS scheme. 

The renormalization constant Z3 is, in the GMS scheme, defined by 

z-^ = 1 + n(//2) (2. 3) 

where n(A;^) is the gluon self-energy appearing in the transverse part of gluon 
propagator 

A correct way of calculating the function Il{k^) is to use the relation 

n(^^) = ^^"'^n,.(A;) (2. 5) 

where n^i,(/c) is the radiative correction tensor which is transverse, Il^j,{k) = 
{k'^g^u — fc//^i/)n(A;^), as implied by the Ward identity k'^Iln^{k) = 0. For 
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the one-loop gluon self-energy diagrams as depicted in Figs.(la)-(lc), after 
computing the Il{k^) by the dimensional regularization procedure and then 
setting = /i^, one can get the regularized form of the Z3 as follows 

-.2 



Z, = 1+ » (4rf)T(l+.)/;d.{^ 



r 

(47^)2 V^—^ ^ -V- ■ -yj^ ^"^ ' [/i2a;(a; - 1)]^ 

-Jlix) + J2(x)] + -y , , , — ^ 



where 



Mx) = 2(n-l)a;(x-l){^-[3n- 4 + 2(1 -O(x-l)] 

n — 2 

J2(x) = -(l-0(ri-l)[4(l + -)x(x-l)(x-2) + (l + -)x 

n n 

(2.8) 

73(0;) = ^n(n - l)a;(a; - 1) (2. 9) 

and rrii is the mass of i-th quark, while, M is an arbitrary mass introduced 
to make the coupling constant g dimensionless in the n-dimensional space. 
The constant Z3 is, in the GMS scheme, defined as 

Zg-^ = 1 + cuifi^) (2. 10) 

where is the self-energy of ghost particle which appears in the ghost 

particle propagator 

i^'^^Q) = (2. 11) 

For the one- loop diagram shown in Fig.(le), the regularized form of the Z3 is 
easily obtained by the dimensional regularization as shown in the following 

Z3 = l + -^(47rM2)T(l + £) C dxm-x){ , , / 

x[^(l + + 3(l-0^] + J(l-O(l-10a:)} (2. 12) 



Now, let us turn to discuss the subtraction of the ghost vertex. The vertex 
is generally represented as 

rf%P,q)^9r''\p^ + Mp,q)] (2. 13) 

where A^(p, q) denotes the higher order corrections. According to the Lorentz 
covariance, it may be written in the form 

= ^(p^?^p•?)PM+^(p^^^p•^)^M (2- 14) 

where A and B are the scalar functions. Suppose only the function A is 
divergent in the limit e ^ 0, while the B is finite as we encountered in the 
one-loop approximation. In this case, we may only subtract the divergent 
part in the A at the point = q^ = /j,^. Alternatively, it is more convenient to 
choose the renormalization point such that p=q and = fx^ which implies 
— {p — qY — 0. This choice is compatible with the energy-momentum 
conservation and yields 

|p=g,p2=^2= Lp^ (2.15) 

where 

L = A{fi^) + B{f?) (2. 16) 

With the subtraction in Eq.(2.15) and the definition 

^ = 1 + Z (2. 17) 

The vertex in Eq.(2.13) may be renormalized as 

^T(P:Q)-zr'rtM<i) (2- 18) 

where 

rt(P> Q) = ar'^iP, + AR,iP, q)] (2. 19) 

which is the renormalized vertex. The latter vertex satisfies the boundary 
condition 

r£b'?) l™^=M^=^r'% (2- 20) 

which is the bare vertex, showing the advantage of the subtraction chosen. 
For the one-loop diagrams depicted in Figs. (If) and (Ig), according to the 
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definition in Eqs.(2.15) and (2.16) and through tedious calculations, we find 
in the n-dimensional space 

L^Li+L2 (2.21) 
where Zi and L2 are given by diagrams in Figs.(lh) and (li) respectively. 

^' = + e) /; d.jpj^-^i(3« - 2 + 3(1 - m + fts 

(2. 22) 

and 

+fts (2. 23) 

where the symbol" fts" represents the terms which are finite in the limit 
£ — > 0. These terms are not necessary to be written explicitly for our purpose 
because they are independent of the renormalization point /i and therefore 
give no contributions to the anomalous dimension. In the approximation of 
order g'^, considering Eq.(2.21), the Zi defined in Eq.(2.17) can be written 
as 

Zi = 1 - Zi - Z2 (2. 24) 



3. Effective Coupling Constant 

The effective couphng constant is determined by the following RGE 

/^^^i?(/^) + 75(/^)^r(a^) = (3. 1) 

which is obtained by differentiating Eq.(2.1) with respect to fi. According to 
the definition in Eq.(2.2), the anomalous dimension 7g(/x) is given by 

7g = J™/^^ In ^9 = 7i - 73 - ^73 (3. 2) 

The anomalous dimensions 73, 73 and 71 in the approximation of order are 
easily calculated from the corresponding renormalization constants. Prom 
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Eqs.(2.6)-(2.9), we obtain 



73 = ii^/^^1^^^3(/^,^) 



n- 



^ :{3e-13 + tE[l + 6^.'+ ' 



xln ^^]} (3.3) 

1 - ./l - Aa} 



where c^i = Prom Eq.(2.12), it follows 

By using the expressions given in Eqs.(2.21)-(2.24), it is easy to find 

d ~ 

There are three points we would like to stress here. (1) The results given in 
Eqs.(3.4) and (3.5) are exactly identical to those obtained in the MS scheme. 
This is because in the one-loop diagrams of the ghost particle self-energy 
and the ghost vertex, only the massless particles are involved. So, the 71 and 
73 given above are scheme- independent. This is why we like to choose the 
ghost particle self-energy and the ghost vertex to define the renormalization 
constant Zg . As for the 73, it is noted that the term related to the J2(a;) 
in Eq.(2.6) which is finite in the limits — > gives no contribution to the 
anomalous dimension owing to its independence of the fermion mass. When 
the quark mass is set to be zero, the result in Eq.(3.3) will be reduced to 
that given in the MS scheme. In this case, the 73 is only given by the terms 
in Eq.(2.6) which contain the e-pole. (2) Each of the quark bare masses in 
Eq.(3.3) has been replaced by its renormalized one which is taken to be a con- 
stant and identified with the pole of the quark propagator. The replacement 
is suitable only in the lowest order approximation. At one-loop level, we con- 
cern the anomalous dimension of the order of ^ therefore, the quark mass 
in Eq.(3.3) is only needed to be given in the lowest order. (3) The renormal- 
ization point may be parametrized by a scale variable A, = /xqA where \Xq 
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is a fixed scale parameter which may be chosen from physical consideration. 
The A can also be chosen to be the scale variable for momenta, p = Xpo. 
With the A introduced, we may write cTj = ^ where q = Noticing this 
relation, on inserting Eqs.(3.3)-(3.5) into Eq.(3.2), the anomalous dimension 
7g will be represented in the following 

(4/1 j 



where 



in which 



F9W = 11 - 3 E[l + ^ + ^/^(A)] (3. 7) 



1 A + JX^ - 4cf 

/,(A) = In V 



- 4cf A - ^A^ - 4cf 



= < 



^4= cot-i , ^, ■ if A < 2q 
^J== coth~^ , ^ if A > 2ci 



(3. 8) 



Substituting Eq.(3.6) into Eq.(3.1) and noticing //^ = A^, the equation 
may be written as 

On integrating the above equation by applying the familiar integration for- 
mulas, we obtain 

an{\) = (3. 10) 

1 + If G(A) 



2 

where aR — ^ and 



47r 



in which 

Nf 



« = E[2c.'-(l + 2cDx.] (3. 12) 



i=l 
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here 

X^ = Vl-4cfln^(l + Vl-4cf) 



4cf - Icos ^ if 2ci>i 



(3. 13) 



^1 - 4cf cosh ^2^, if2Q<i 

& = Ec' (3.14) 



and 



here 



=1 



^W = E(1 + ^)X^^W (3- 15) 



r,,{\) = VA2-4cfln— (A + VA2-4c|) 



4c| - A^cos-^ 2^, if A < 2c 



^A2 -4cf cosh-^ ^, if A > 2q 



(3. 16) 



If we set all the quark masses to be equal and choose /xq = rnji, Eqs.(3.7) 
and (3.11) will be respectively reduced to 



and 



F,(A) = 11 - lNf[l + 4 + In ^±^^] (3. 17) 



G(A) = 11 In A - ^7V^[2 + v^Stt - ^ + (1 + ;|)^ry(A)] (3. 18) 



where 



77(A) = ^A2^1n^(A + VA2^) 

f-V^^^cos-if, ifA<2 
I ^/F^cosh-^ |, ifA>2 



(3. 19) 



In the above formulas, the A may be defined by A = |g^/m|.|2. When A = 1, 
due to G(1)=0, Eq.(3.10) becomes ^^(1) = ct^j- This gives the result on the 
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fermion mass shell. When A — > cxo, i.e. in the large momentum limit, it is 
easy to sec from Eqs.(3.18) and (3.19) that G(A) ^ (11 - |iV/)lnA. Thus, 
Eq.(3.10) is reduced to 

an{X) = ^^r-^ ^ (3. 20) 

l + ^(33-2iV^)lnA 

This just is the result obtained previously in the MS scheme^'^. The behavior 
of the effective coupling constants arc represented in Figs. (2) and (3). The 
coupling constants in Fig. (2) are given by taking all the quark masses to 
be equal and Nf = 4. The solid curve represents the coupling constant 
evaluated by using Eqs.(3.10), (3.18) and (3.19). Fig. (3) shows some coupling 
constants for different numbers of flavor and unequal masses. In the flgure, 
the constituent quark masses are taken to be m„ = rrid = 0.3GeV,ms — 
0A5GeV,mc = 1.5GeV,mi) = 5.0GeV and rrit = 175GeV, and the scale /iq 
is taken to be equal to the A^/-th quark mass. From Figs. (2) and (3), one 
can see that all the q;_r(A) decrease with the increase of A and tend to zero 
when A ^ oo, exhibiting the well-known asymptotically free property. In the 
region near A = 1, there is a maximum for each curve. The height and the 
position of the maximum weakly depend on the number of flavor(see Fig. (3)). 
In comparison with the result given in the MS scheme, as we see from Fig. (2), 
the quark mass gives a considerable improvement on the effective interaction, 
particularly, in the region near the heavier meson threshold. However, as 
shown in Figs. (2) and (3). When A — > 0, all the Q!i?(A) drastically fall down 
to zero from their maxima. This unreasonable result indicates that the QCD 
perturbation theory is inapplicable in the very small momentum domain. 



4. The Ward Identity for Quark-Gluon Vertex 

The aim of this section is to sketch the derivation of the Ward identity 
satisfled by the quark-gluon vertex and to show how we should do for the 

subtraction of the quark self-energy. The Ward identity is not difficult to 
derive from the QCD generating functional. Firstly, we write the Ward 
identity obeyed by the quark-gluon three point Green function^^ 

< 0+\T[^{x)^(y)Al(z)]\0-> 
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= tig{< 0+ I T[i,{x)i;{y)C\y)c\z)] | 0" > T" 
-T^ < 0+ I T[i;{x)^{y)C\x)d'' {z)] \ 0" >} (4. 1) 

where denote the generators of SU(3) group, ■0(x),74^(x) and C"(x) 
represent the field operators for quark, gluon and ghost particle respec- 
tively. The Green's functions in Eq.(4.1) have the following irreducible 
decompositions^^'^^ 

< 0+ I T[ip{x)^{y)Al{z)] I 0- > 

= j d'x'd'y'dh'SFix - x')T'^{x\ y\ z')SF{y' - y)D':^{z' - z) (4. 2) 

< 0+ I T[i,{xmy))C\y)t\z)] \ 0" > 

= j d^x'd^z'Spix - x')7^"' {x', y, z') A""' {z - z') (4. 3) 

where 



ha! I t I 

7 K^^V^z . 



J d'y'd'u'T''^'{x', y', u', z')Sp{y' - y)^'>'\u' - y) (4. 4) 
and 

< 0+ I T[ip{x)1p{y)C\x)d\z)] I 0" > 

= / d'y'd'z'-i^'^'ix, y\ z')SF{y' - A"''^(^' - z) (4. 5) 

where 

7^"'(x, y, z') = j d^x'd^u'Spix - x')r''''''{x', y, u', z')A'''\u' - x) (4. 6) 

In the above, Sf{x — x'), D'^'^{x — x') and A"^(x — x') are the quark, gluon 
and ghost particle propagators respectively, r°^(a;, y, z) denotes the three- line 
quark-gluon proper vertex and r''^(a;, -u) designates the four-line quark- 
ghost vertex. Upon substituting Eqs.(4.2), (4.3) and (4.5) in Eq.(4.1) and 
transforming Eq.(4.1) to the momentum space, the Ward identity will be 
written as 

SF{p)T'''{p,q,k)SF{q)k^DfM 
= ^9b\p, k)SFiq) - 5f (p)f (p, q, k)]A^\k) (4. 7) 
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where 

^\p,q,k) = T^^''^ip,q,k) (4. 8) 

i>ip,q,k) = ^"^ip,q,k)T'^ (4. 9) 

Employing the expressions denoted in Eqs.(2.4) and (2.11) and operating on 
the both sides of Eq.(4.7) by Sp^{p) from the left and by Sp^{q) from the 
right, we arrive at 

k^ri{p,q,k) = g[l + Lo{e)]-'[S^\pmP,q,k) 

-riP,q,k)S^\q)] (4. 10) 

This just is the Ward identity satisfied by the quark-gluon vertex^^. Consider- 
ing the energy- momentum conservation and introducing new vertex functions 
A^(p, q) and x"(p, q) which are defined by 

T;ip, q, k) = i27iyd\p - q - k)gA;ip, q) (4. 11) 

nP,q,k)^(27ry5\p-q-k)[l + u;(k')]x''(p,q) (4. 12) 

r(p,q,k)^(2ny5\p-q-k)[l + u;(e)]r(p,q) (4. 13) 

the identity in Eq.(4.10) can be represented in the form 

(p - qTKip^ q) = s-p'{p)x''{p, q) - rip, q)sAq) (4- 14) 

In the lowest order approximation, the above identity is clearly satisfied as 
long as we notice 

Ai°^'^(p,9) = ^7/.T« (4. 15) 

X(°)«(P,?) = X^°^"(p,?)=^r« (4. 16) 

The latter result may easily be derived from the definitions given in Eqs.(4.4), 
(4.6), (4.8), (4.9), (4.12) and (4.13) by using the lowest order expression 

V^'ix, y, z, u) ^ iS%{x - y)-'^l,{z - u)-' (4. 17) 

Now, we are interested in the one-loop approximation of order g^. In 
this approximation, the quark-gluon vertex denoted by A^^^"(p, g) is con- 
tributed from the two diagrams in Figs.(lh) and (li) whose expressions are 
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well known. The quark-ghost vertex functions q) and x"(p, q) evaluated 
from Figs.(lj) and (Ik) are shown in the following. 

X^'^''{p,q)^iT''I{p,q) (4. 18) 

where 

I{p,q) = -^g'l ^^S^{l)rD,.{q - l){l - pT A{1 - p) (4. 19) 

and 

X^'^^{p,q) = iT^7{p,q) (4. 20) 

where 

Hp,<i) = -la'J^.YSFmi-qr 

xD^,{l-p)A{l-q) (4. 21) 

It is clear that the above functions are logarithmically divergent. Thus, up 
to the order g^, we can write 

A'i^%p,q) = tT'^b, + A'i\p,q)] (4. 22) 

where we have set A^(p, q) — iT°'A.^^\p, q), 

X»(p,g)=^T''[l + /(p,g)] (4. 23) 

and 

r(p,g)=ir«[l+7(p,g)] (4. 24) 

Upon differentiating the both sides of Eq.(4.14) with respect to p^, set- 
ting q=p and substituting the expression of the inverse of quark propagator 

Sp\p) = i)-m-i:{p) (4. 25) 

in the order of g"^, we get 

A,(p,p) = (4. 26) 
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where 



I ( If (1 ] 

I ,.p(/-m) (4.27) 



^p^^ 

It is emphasized that at one-loop level, the both sides of Eq.(4.26) are of the 
order of g"^. In the derivation of Eq.(4.26) from Eq.(4.14), the terms higher 
than the order have been neglected. The identity in Eq.(4.26) formally is 

the same as we met in QED. The vertex A^(p,p) may be expressed in the 
form 

A^(p,p)=L7^ + A^(p) (4. 28) 

where L is a divergent constant defined by 

L = A^(p,p)|j,=^ (4. 29) 

and A^(p) is the finite part of A^(p,p) satisfying the boundary condition 

A^(p) |^=M= (4. 30) 

Integrating the identity in Eq.(4.26) over the momentum and considering 
the expression in Eq.(4.28), we have 

E(p) = A + (i>-^)[5-C(p2)] (4. 31) 

where A and B are the divergent constants depending on the renormalization 
point /X which are defined as 

A = E(/x) (4. 32) 

B^-L (4. 33) 

and C(p^) is a finite function defined by 



with boundary condition 



rdp'^Al{p)^{^-pi)C{p') (4. 34) 

C(p2) = (4. 35) 
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Clearly, the expression in Eq.(4.31) gives a subtraction version for the fermion 
self-energy which is required by the Ward identity and correct at least in the 
approximation of the order g^. With this subtraction, the quark propagator 
will be renormalized as follows 

Sf{p) = ^ ^ (4. 36) 

where Z2, nii^ and 'Er{p) are the quark propagator renormalization constant, 
the renormalized quark mass and the finite correction of the self-energy re- 
spectively. The Z2 and the vhr are respectively defined as 

Z-^ = 1-5 (4. 37) 

and 

ruR = Z-^m (4. 38) 

where is the quark mass renormalization constant which is defined by the 
following expression 

= 1 + Z2[Am-^ + (1 - /xm-^)5] (4. 39) 



5. Effective Quark Mass 

Taking the derivative of Eq. (4.38) with respect to n and noticing = A^, 
we get a RGE for the renormalized quark mass as follows 

\^mR{\) + 7m(A)mi^(A) = (5. 1) 

OA 

where ^ 

7m(A) = lim/x— lnZ„ (5.2) 

is the mass anomalous dimension. Let us concentrate our attention on the 
one-loop approximation. The fermion self-energy of the one-loop diagram 
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shown in Fig. (11) is of the following regularized form in the n-dimensional 



spacc'^^'^^ 



x[2(l -s){l-x) i)- 2(2 - e)m + {I - ^){m - 2xp)\ 
-2(1 - 0(1 -x)^^^!^} (5- 3) 



A(p2) 



where 



A(p^) = p^x{x — 1) + rn^x (5. 4) 

Substituting the above expression in Eq.(4.31), one may find 

X [2^l[l + (C - 2)x - £(1 - x)] - (3 + ^ - 2e)m] 
-2(1-0(1-^)^^^^} (5.5) 



and 



K4.M^)T(l + s)/;d.{-l 



127r2' ' ^ ' Jo '£A(/i2)e 

2u2 

x[2(l - e){l -x)- 2(1 - Ox] + -^^A^ - 1) 

777 

X [2(1 -e){x-l) + 5(1 - Ox + -(3 + e - 2e)] 

A* 

(l-0(l + £)(x-l)V} (5.6) 



A(/.2) 



2+£' 



where 



A(/x2) - x[{x - 1)/^^ + m^] (5. 7) 

By making use of the renormalization constants defined in Eqs.(4.39) and 
(4.37) and the constants given in Eqs.(5.5) and (5.6), in the order of g^, it is 
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not difficult to derive the anomalous dimension defined in Eq.(5.2) for i-th 
quark as shown in the following 



7^^ (A) 



TT 



(5. 8) 



where 



|^A + 2[3 + 2e-3(l + Oy 



A2 ^ 

-3(1 + 0y 



Ci + X 



A2 



In 



c? - A2 



(5. 9) 



in which a — The //q may be chosen as stated in Sect. 3. As seen from 
Eq.(5.9), the anomalous dimension and hence the effective mass are gauge- 
dependent at one-loop level. In the following, we arc only interested in the 
result given in the Landau gauge which was regarded as preferred gauge in 
the previous literature^''. In this gauge, 



X In 



1 - 



4A 6c? _ Ci 
A + Q A2 ^ A^ 
A^ 



(5. 10) 



Substituting Eqs.(5.8) and (5.9) in Eq.(5.1) and solving the equation, the 
effective mass is found to be 



mg(A) 



where 



^«(A) = -/'^MA)F«(A) 
n Ji A 



(5. 11) 



(5. 12) 



here the q;r(A) was shown in Eq.(3.10). If we take aR{X) ~ and work in 
the Landau gauge, we get 

T ,2c? 3c2 



^«(A) = v^2Q(i-i) + (^-:^ + i] 

X In 1 5" 

cf 



- {24 - 3c^ + 1) In 



1 - 



} (5. 13) 
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If we take q = 1, namely, set all masses to be equal and choose /iq = rrf^, 
Eq.(5.13) will be reduced to 



al (1 - A) 




1 + A 
A 



] In 1 - A^ } 



(5. 14) 



TT A 



where 



In 1 - A^ 



f2[ln(l + A) -tanh-^A], ifA < 1; 
|2[ln(l + A) -coth-^A], if A > 1 



(5. 15) 



The behavior of effective masses is illustrated in Fig. (4). The results in 
Fig. (4) are given by taking Nj = 4 and all the quark masses to be equal. 
The solid curve and the dashed one represent the effective masses calculated 
respectively by taking the coupling constant to be the running one and a fixed 
value aR = 0.2. The common feature of these effective masses is as follows. 
For each curve, there is a maximum at A ~ 1.54. When A tends to infinity, 
the masses fall down to zero, exhibiting the well-known asymptotically free 
behavior. For the A lying in the region [0,1), the mass is less than the 
maximum and behaves almost as a constant. However, in the region of A 
near zero, as pointed out in Sect. 3, the QCD perturbativc results are no 
longer valid. 

6. Comments and Conclusions 

In this paper, the QCD one-loop renormalization has been restudied in 
the GMS scheme. It was shown that the GMS scheme allows us not only to 
consider the mass effect on the renormalization, but also to directly relate the 
renormalization scale /i to the momentum p. The effective coupling constant 
and the effective quark masses obtained in this scheme get a noteworthy im- 
provement near the heavy quark threshold in comparison with those given 
previously in the MS scheme. The effective coupling constant and the ef- 
fective quark mass presented in the MS scheme now appear as approximate 
results given in the large momentum limit. It is noted here that even in 
this limit, the effective coupling constant and the effective mass can only 
have unique forms. It is impossible to result in a difference between the MS 
and MS schemes^'®'^^ because an effective quantity should be the solution 
of its RGE whose form is uniquely determined by the anomalous dimension 
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and the boundary condition. As emphasized before, the correct procedure 
of evaluating the anomalous dimension is the use of the regularized renor- 
malization constant and the limit £ — > should be taken after differentiation 
with respect to the renormalization point. By this procedure, the factor 
(47rAf)'^r(l + e) appearing in the n-dimensional Feynman integrals can only 
come to unity. It is not possible to yield the unphysical terms In 47r — 7 in a 
renormalized quantity. 

In this paper, the effective coupling constant given in the one-loop ap- 
proximation is calculated by employing the renormalization constants given 
by the subtractions of the gluon and ghost particle self-energies and the ghost 
vertex correction. It is emphasized that the renormalization constant Z2 is 
obtained on the basis of the Ward identity obeyed by the gluon self-energy 
and therefore is faithful. The renormalization constants Zi and Z3 derived 
respectively from the subtractions of the ghost vertex correction and the 
ghost particle self-energy, as mentioned in Sect. 4, give the anomalous dimen- 
sions which are scheme-independent due to that there is no quark mass to 
appear in the vertex and the self-energy. Therefore, the anomalous dimen- 
sion 7g computed from these renormalization constants is definite, no any 
uncertainty. Particularly, the correctness of the anomalous dimension shown 
in Eqs.(3.6)-(3.8) is confirmed by the previous result presented in Refs.(9) 
and (16) where the 7^ was determined by the subtraction of the quark and 
gluon self-energies and the quark-gluon vertex at a space- like point. The 
expression of the /3- function given in the space-like point is 

A3(A2 + 4cf)^ (A2 + 4c|)i-A 

where we have set /i — /IqX and q — —. The above expression may be 
directly written out from Eqs.(3.6)- (3.8) by the transformation c? — > — c?, 
corresponding to /x^ — > — /i^. The effective coupling constant given by the f3- 
function in Eq.(6.1) still exhibits the property of asymptotic freedom when 
A — > cxd; but, it has a Landau singularity in the region of large distance, 
different from the result given in the GMS scheme. In Fig. (2), the Landau 
singularity occurs at A 0.4. 
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It is stressed again that the effective quark masses shown in Sect. 5 are 
obtained based on the subtraction written in Eq.(4.31) which is derived from 
the Ward identity respected by the quark one- loop self-energy. These effective 
masses are different from those given in the previous works^'^''. For example, 
in Ref.(3), the mass anomalous dimensions derived in the Landau gauge is 
of the form 

which is manifestly different from that formulated in Eqs.(5.8) and (5.10). 

The result shown above , actually, was obtained from another subtraction. 
If we start from the expression for the quark self-energy^^ 

E(p) = A{p^)P + B{p^)m (6. 3) 

and subtract at the space-like point = —fj.^, the quark propagator and 
mass renormalization constants will be defined by 

Z,' = 1 - A{-^') (6. 4) 

and 

Z-' = Z,[l + Bi-fi')] (6. 5) 

In the approximation of order for the one-loop diagram, the constant 
defined in Eq.(6.5) will give rise to the result in Eq.(6.2). Since the sub- 
traction above is not compatible with the Ward identity, the effective quark 
masses obtained in this subtraction can not be viewed as reasonable results. 
The correct result of the effective quark mass given by the subtraction at 
space-like may be obtained from Eqs.(5.11)-(5.15) through the transforma- 
tion A — > iX. In this case, as we see, the effective mass becomes complex. If 
we require the effective mass to be real, the subtraction at space-like point 
should be excluded. However, this does not mean that the space-like momen- 
tum subtraction is useless. The subtractions at time-like point and at space- 
like point probably suit to different processes in which the interactions are of 
different behaviors. For example, for the t- channel fermion-antifermion scat- 
tering, the momentum on the intermediate boson line is space-like; while, for 
the s-channel scattering, the corresponding momentum is time-like. It seems 
that the boson self-energy should be subtracted at the space-like point for 
the former process and at time-like point for the latter process. This problem 
is, we think, worthy to pursue in future investigations. 
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8 FIGURE CAPTIONS 



Fig.(l) The one-loop diagrams. 

Fig. (2) The one-loop effective coupling constants given by taking = 
0.2, Nf = 4 and equal mass of quarks. The solid curve represents the result 
obtained at time-like subtraction point. The dashed curve represents the 
one given at space-hke subtraction point. The dotted curve shows the MS 
scheme result. 

Fig. (3) The one-loop effective coupling constants evaluated for the dif- 
ferent numbers of flavors Nj = 1,3,6 and for unequal quark masses. The 
masses of quarks of different flavors are taken to be the values as written in 
Sect. (3). 

Fig. (4) The one-loop effective quark masses obtained by taking Nf — A 
and equal mass. The solid curve represents the effective mass obtained by 
using the running coupling constant. The dashed curve represents the one 
given by taking the coupling constant to be a constant a^j = 
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Fig. 3 
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